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MEMORANDUM 


To:  All  Authorized  Recipients  of  the  DNA  Reaction  Rate 

Handbook  (DNA  1948H) 

From:  The  Editors 

Subject:  Revision  Number  2 


Enclosed  herewith  you  will  find  a  copy  of  Revision  Number  2  to  the 
Handbook.  It  comprises  the  newly  written  Chapter  22.  Related 
revisions  to  appropriate  pages  of  the  introductory  matter  and  the 
Appendices  will  be  included  in  a  later  revision. 

You  should  immediately  substitute  the  enclosed  items  in^o  your  copy 
of  the  Handbook,  discarding  old  pages  22-1  and  22-2. 

You  should  also  enter  on  page  iii  in  front  of  your  Handbook  the  follow¬ 
ing  information:  Revision  No.  2;  Date  of  Issue— May  1974;  Date  of 
Receipt  — whatever  day  you  receive  this;  and  sign  your  name  in  the 
last  column. 


Revision  Number  3  has  already  been  issued.  Revision  Number  4  is 
expected  to  be  issued  during  the  late  summer  or  early  fall  of  1974. 
Thank  you  for  your  cooperation. 


It  will  be  noted  that,  as  these  revisions  continue  to  be  issued,  some 
of  the  Handbook  appendices  will  tend  to  become  outdated.  From  time 
to  time  revisions  to  the  appendices  will  be  distributed,  in  order  to 
correct  this  situation. 
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CHAPTER  22 


22.  APPLICATIONS  OF  COMPUTER  SOLUTIONS  TO 
ATMOSPHERIC  DEIONIZATION  PROCESSES 

A.W.  All,  Naval  Research  Laboratory 
W.S.  Knapp,  General  Electric— TEMPO 
F.E.  Niles,  Ballistic  Research  Laboratories 

(Latest  Revision  5  March  1974) 


22.1  INTRODUCTION 

Many  natural  and  man-made  events  entail  disturbances  of  the 
upper  atmosphere.  Of  particular  interest  is  the  sequence  of  deioni¬ 
zation  phenomena  following  a  nuclear  burst.  The  radiation  and 
debris  from  a  nuclear  burst  interact  with  the  ambient  atmosphere  in 
the  three  environmental  regions  listed  in  Chapter  1,  viz.  ,  (1)  the 
fireball  region,  (2)  the  fireball  halo  or  the  UV  fireball,  and  (3)  the 
■~ool  disturbed  region  beyond.  This  interaction  results  in  a  disturbed 
atmosphere  with  varying  degrees  of  ionization  and  additional  species 
which  are  not  present  or  are  very  minor  under  normal  conditions. 
These  disturbed  regions  can  be  described  mathematically  with  an 
appropriate  computer  code  designed  for  a  specific  function.  In  gen¬ 
eral  one  is  interested  in  the  time  histories  of  the  species  in  the  dis¬ 
turbed  regions  and  the  time  histories  of  the  electron  temperature, 
the  vibrational  temperature,  and  the  heavy-particle  (neutral)  tempera¬ 
ture.  Thus  rate  equations  are  written  for  each  species  and  the  de¬ 
sired  temperatures.  Then  these  equations  are  solved  using  some 
type  of  numerical  integration  technique  to  obtain  the  desired  quantities 
as  a  function  of  cime.  Codes  with  many  species  are  known  as  multi- 
s  pecies  codes  . 

However,  to  describe  the  deionization  of  a  disturbed  region  by 
an  appropriate  code  one  must  know  the  source  ol  the  disturbance, 
its  duration,  the  atmospheric  species  densities,  and  basic  physical 
and  chemical  processes  involved.  (See  Chapter  6  for  a  discussion  of 
the  chemical  kinetics  of  the  disturbed  atmosphere.) 

On  the  other  hand  the  complexity  of  a  code,  e.  g.  ,  the  number  of 
species,  temperatures,  reactions,  and  other  desired  quantities  de¬ 
pends  on: 

(a)  The  purpose  for  which  the  code  is  designed,  i.e.  ,  whether 
the  code  is  to  be  used  for  radar  propagation,  communication, 
or  optical  systems  studies,  or  for  some  or  all  of  these  pur¬ 
poses  simultaneously. 
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(b)  Whether  the  code  is  to  be  used  in  conjunction  with  hydro- 
dynamic  codes. 

(c)  Computational  time  requirements  and  the  size  of  the  com¬ 
puter,  since  the  time  required  to  perform  detailed  deioniza¬ 
tion  calculations  could  be  prohibitive  where  solutions  are 
needed  at  thousands  of  points  in  space  and  time. 

No  computer  code  can  include  all  of  the  possible  reactions  which 
may  occur,  especially  when  each  energetic  state  of  a  constituent 
should  be  considered  as  a  separate  species  with  its  own  reactions. 
One  hopes,  however,  that  the  principal  species  and  reactions  aie 
identified  and  included.  In  this  chapter  several  multispecies  codes 
for  the  D,  E,  and  F  regions  are  discussed.  The  dominant  reactions 
and  the  important  species  are  pointed  out  for  each  region. 

Numerical  methods  for  solving  the  rate  equations  are  reviewed. 
Simplified  codes  with  physical  and  mathematical  approximations  and 
the  lumped- parameter  method  are  also  discussed. 


22.2  NUMERICAL  METHODS  FOR  SOLVING 
THE  RATE  EQUATIONS 

The  rate  equations  for  n  species  are  a  set  of  n  first-order, 
coupled,  and  generally  nonlinear  differential  equations.  Given  a  set 
of  chemical  reactions  involving  the  n  species  and  specifications 
for  the  concentrations,  nQ  at  time  tg  for  the  i'th  species,  the 
requirement  becomes  one  of  generating  a  set  of  n  differential  equa¬ 
tions  and  solving  the  set  numerically  for  time  greater  than  t^  . 

The  differential  equations  may  be  written  as: 


(22-1) 


where  fj  is  the  rate  of  formation  of  the  i'th  specie9  and  (t^)*!  is  the 
removal  frequency  of  the  same  9pecles, 


If  fj  and  t-  are  constants,  Equation  (22-1)  can  be  solved  analyti¬ 
cally  to  yield: 

n.(t)  *  i  (nQ  .  -  f.t.)  exp  f-(t  -  t^/tj  .  (22-2) 
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Equation  (ZZ-2)  reveals  that  1/t^  describes  how  quickly  the  concen¬ 
tration  of  the  i'th  species  reaches  its  equilibrium  value  of  fj  tj  with  a 
characteristic  time  of  t^  . 

Unfortunately  fj  and  t^  are  normally  functions  of  time  and  the 
differential  equations  most  often  encountered  have  time  constants 
for  the  different  species  which  vary  over  many  orders  of  magnitude. 
Often  the  shorter  time  constants  are  much  shorter  than  the  time 
step  desired  for  performing  the  numerical  integration.  For  the 
classical  methods  of  numerical  integration,  the  solutions  become 
unstable  if  for  the  largest  (t^)-^  ,  denoted  as  R,  the  product  of  R 
and  the  time  step  of  integration,  h,  is  greater  than  some  small 
quantity  which  is  typically  near  one,  i.e.  ,  Rh  >  1.  The  equations 
for  which  (t^)"^h  is  £  1  are  said  to  be  numerically  stiff.  When  this 
is  the  case,  the  problem  of  integration  becomes  more  complex  and 
nonclassical  methods  (which  are  free  of  the  constraints  that 
(tj)~^h  <  1)  must  be  adopted  in  order  to  obtain  solutions  -within  a 
reasonable  amount  of  computer  time.  The  nonclassical  methods 
must  confront  two  major  restraints,  viz.  ,  accuracy  and  stability. 

The  nonclassical  methods  fall  into  two  classes:  those  which  treat 
all  the  differential  equations  in  the  same  manner  and  those  which 
recognize  the  stiff  equations  and  integrate  them  by  a  stable  method 
with  the  rest  o:  the  equations  being  integrated  by  a  more  accurate, 
but  less  stable,  classical  method.  We  shall  refer  to  the  latter 
class  as  "hyarid  methods.” 

In  all  hybrid  methods  two  important  considerations  mu3t  be  ad¬ 
dressed.  The  first  is  the  establishment  of  a  criterion  to  determine 
by  whicu  method  specific  equations  are  to  be  integrated.  The  second 
is  hc»v.  the  chemical  balance  and  charge  balance  are  to  be  maintained 
and  how  much  imbalance  will  be  allowed.  Hybrid  methods  are  rea¬ 
sonably  fast;  however,  the  accumulated  errors  may  become  large 
ait«  r  many  time  steps  so  that  care  must  be  exercised  in  using  them, 
ih.neshea  (Reference  22-1)  has  developed  and  utilized  a  hybrid 
ethod  in  his  codes. 

The  main  ideas  behind  the  classical  methods  along  with  a  survey 
and  discussion  of  the  literature  have  been  given  by  Benyon  (Reference 
22-2).  He  has  given  comparisons  for  predictor,  predictor  corrector, 
and  Runge-Kutta  methods  and  has  listed  75  refe rence s .  Also  he  has 
provided  a  few  comments  regarding  nonclassical  methods.  Several 
books  have  been  written  on  methods  of  solving  differential  equations. 
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A  few  of  these  are  listed  among  the  references  at  the  end  of  this 
chapter  (References  22-3  through  22-21). 

Among  the  nonclassical  methods  of  principal  interest  for  atmo¬ 
spheric  deionization  calculations  are:  the  hybrid  method  of  Keneshea 
(Reference  22-1),  the  algorithm  of  Treanor  (Reference  22-22),  the 
explicit  method  of  Fowler  and  Warten  (Reference  22-23),  the  selected 
order,  multistep  predictor -corrector  method  of  Gear  (Reference 
22-241,  the  matrix  method  of  Kregel  (Reference  22-25),  and  the 
CHEMEQ  Stiff  Equation  Integrator  by  Young  and  Boris  (Reference 
22-26). 

Young  and  Boris  (Reference  22-26)  have  devised  their  method  for 
the  purpose  of  extremely  fast  numerical  solution  in  the  chemistry 
rate  equations  of  the  deionization  type  while  maintaining  a  moderate 
level  of  accuracy,  and  for  coupling  rate  equations  to  the  general 
flows  arising  in  disturbed  atmospheric  phenomena. 

In  selecting  the  method  to  be  used,  one  obviously  wants  that  which 
will  require  the  least  amount  of  computer  time  while  maintaining  an 
acceptable  accuracy,  A  direct  comparison  of  methods  is  complicated 
by  the  fact  that  they  may  not  all  be  optimized  in  the  coding  for  the 
computer  being  used  for  the  comparison  study.  A  poor  method  for  a 
given  application  may  run  faster  than  a  better  method  when  the  poor 
method  is  coded  using  optimum  coding  practices.  Examples  of  non- 
optinmm  coding  practices  include:  (1)  using  doubly  subscripted 
variables  when  a  singly  subscripted  variable  could  be  used;  and  (2) 
passing  information  in  and  out  of  subroutines  of  functions  using  the 
linkage  rather  than  the  use  of  COMMON  when  COMMON  could  be 
used. 

In  addition  to  considering  which  method  of  numerical  integration 
is  best  and  the  use  of  optimized  coding  practices,  one  must  consider: 
(11  how  to  compute  or  estimate  the  error  per  time  step;  (2)  criteria 
for  changing  step  size;  (3)  how  much  the  step  size  should  be  changed 
when  it  is  changed;  (4)  how  to  conserve  chemical  balance;  (5i  how  to 
conserve  charge  balance;  and  other  factors.  In  many  cases  these 
considerations  can  be  addressed  only  by  trial  and  error. 

All  methods  which  do  not  conserve  chemical  balance  and  charge 
balance  automatically  may  be  made  to  do  so  by  at  least  two  techniques. 
First,  balance  can  be  invoked  by  summing  the  various  concentrations 
to  find  the  errors  and  then  distributing  the  errors  throughout  the 
densities  in  such  a  way  that  balance  is  conserved.  The  major  fault 
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with  this  technique  is  that  one  is  incorporating  a  portion  of  the  errors 
into  concentrations  from  which  the  errors  may  not  have  arisen.  A 
second,  better  method,  but  one  which  is  not  always  practical,  is  to 
take  the  concentrations  generated  at  each  step  and  use  these  as  input 
into  a  final  step  which  does  conserve  balance  automatically. 

The  success  of  multispecies  codes  for  atmospheric  deionization 
depends,  in  addition  to  an  appropriate  description  of  the  physics 
and  chemistry  of  the  ionized  atmosphere,  on  having  fast,  accurate 
methods  of  numerical  integration. 

22.3  MULTISPECIES  CODES 

22.3. 1  The  Disturbed  E  and  F  Regions  Multispecies  Codes 

To  describe  the  deionization  in  these  disturbed  regions,  it  is  in¬ 
structive  to  know  that  the  dominant  ambient  constituents  of  the  E  and 
F  regions  are  O,  N£,  and  C^.  The  disturbing  source,  whether  it  is 
radiation  or  energetic  particles,  will  create  electrons  and  O4 ,  N^, 

O^,  and  N+  as  the  dominant  ions.  The  last  ion  arises  either  from 
dissociative  ionization  of  or  by  the  ionization  of  the  products  of 
the  dissociative  recombination  of  N£  .  This  last  ionization  process 
is  generally  the  dominant  source  of  N*  during  the  presence  of  the 
disturbing  radiation.  Minor  species,  notably  NO,  are  also  ionized. 
Therefore,  in  principle  one  may  write  rate  equations  for  all  the 
above  species,  include  physical  and  chemical  processes,  use  the 
appropriate  rate  coefficients  (see  Chapter  24),  assume  a  tempera¬ 
ture,  and  thus  have  a  deionization  code.  The  question  to  ask  at  the 
outset,  however,  is  whether  this  code  would  give  the  right  electron 
density.  To  answer  this  question,  let  us  examine  the  basic  deioniza¬ 
tion  processes  which  occur.  The  following  dissociative  recombinations 

+  e  -  N  +  N  (22-3) 

Oi  i  e  “*  O  •»  O  (22-4) 

NO+  +  e  -  N  l  O  (22-5) 

are  the  most  important  processes  which  remove  electrons.  The  first 
two  reactions  occur  at  early  stages  of  the  deionization  while  reaction 
(22-5)  becomes  important  at  late  times.  These  reactions  depend  on 
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the  electron  temperature  {Reference  22-27)  as  well  as  the  ion-kinetic 
temperature  (Reference  22-28)(actually  the  vibrational  temperature). 
The  next  important  reactions  are  the  processes  which  remove  an 
atomic  ion  and  replace  it  by  a  molecular  ion,  i.e.  ,  ion-molecule  re¬ 
arrangement  or  charge-exchange  collisions  such  as: 

o‘  *  K'2  -  NO"  r  N  (22-6) 

N"  *  O,  -  NO  *  O 

2 

-  Oz  4  N  (22-7) 

o'  ■  O  -  O*  *  O  (22-8) 

2  2 

N  ♦  NO  **  No"  4  N  (22-9) 

Reaction  (22-6)  is  widely  studied  and  it  depends  both  on  the  nitrogen- 
vibrational  temperature  (Reference  22-29)  and  the  ion-kinetic  tem¬ 
perature  (Reference  22-30),  In  a  completely  dissociated  and  ionized 
region,  on  the  other  hand,  where  only  atomic  ions  are  present,  col- 
lisional  radiative  recombination  is  the  only  available  deionization 
route.  The  processes  discussed  above  indicate  the  importance  of  the 
electron,  the  nitrogen- vibrational,  and  the  heavy-particle  tempera¬ 
tures.  Reaction  (22-9),  for  example,  depends  on  the  formation  of 
NO,  which  in  turn  depends  on  the  heavy-particle  temperature  through 
the  reaction  of  ground-state  atomic  nitrogen  and  molecular  oxygen 
(see  Chapter  24).  However,  the  formation  of  NO  proceeds  faster  if 
the  reactant  atomic  nitrogen  was  in  the  (^D)  excited  state  (see  Chapter 
24).  This  illustrates  the  importance  of  metastable  states  in  the  de¬ 
ionization  process. 

In  addition  to  the  importance  of  N(^D)  in  the  formation  of  NO,  it  is 
also  a  source  of  electron  heating  by  superelastic  collisions.  The 
same  can  be  said  for  0(  *  D)  which  is  additionally  quenched  by 
(Reference  22-31)  with  almost  all  its  internal  energy  being  converted 
(Reference  22-32)  into  kinetic  energy  and  not  into  vibrational  energy 
of  N2  as  was  speculated  earlier  (Reference  22-33).  However,  a 
small  amount  of  the  0(*D)  energy  may  still  end  up  in  the  vibrational 
mode  of  (Reference  22-34). 
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These  metastable  states  will  heat  the  electrons  via  superelastic 
collisions  and  consequently  raise  the  nitrogen-vibrational  tempera¬ 
ture.  This  leads  to  a  higher  rate  for  the  formation  of  NO+  according 
to  reaction  (22-6)  and  hence  faster  removal  (Reference  22-35)  of 
electrons  (Figure  22-1).  The  following  metastable  states  play  a  simi¬ 
lar  role:  0(1S),  N(2P),  Q4(2D),  04(2P),  N*  ^D),  N^S),  O^a1^), 
and  C>2(bl2).  However,  one  may  neglect  0(*S),  N(2P),  0+(^P), 

N*  (*D),  and  N4  (*S)  if  the  purpose  of  the  code  is  the  calculation  of 
electron  density.  On  the  other  hand,  they  cannot  be  ignored  in  opti¬ 
cal  and  infrared  calculations.  O1  (^D)  plays  an  important  role  since 
it  resonantly  charge-exchanges  (Reference  22-36)  with  : 

o'(2D)  !  N2  -  N2(A2n)  +  O  ,  (22-10) 

and  thereby  changes  an  atomic  ion  into  a  molecular  ion  with  a  rate 
coefficient  of  ~(1  —  3)  x  10"^  cm^  sec'*  (Reference  22-37). 

To  answer  the  question  raised  earlier  on  the  accuracy  of  the  elec¬ 
tron  density  calculations,  it  may  be  concluded  from  the  processes 
discussed  above  that  one  must  include  important  metastable  states  or 
at  least  account  for  them  in  a  realistic  manner.  In  addition,  the  ap¬ 
propriate  temperatures  must  also  be  alculated  as  functions  of  time. 

22.3.2  Examples  af  Deionization  Coleulotions  in 

Disturbed  E  and  F  Regions 

The  Naval  Research  Laboratory  (NRL)  Master  Deposition  and  Late 
Time  Chemistry  Code  (References  22-38  through  22-40)  (hereinafter 
referred  to  as  the  NRL  Master  Code)  and  The  Simple  Code  (Refer¬ 
ence  22-40)  describe  the  deionization  processes  in  disturbed  E  and  F 
regions  of  the  ionosphere.  There  are  other  deionization  codes  (Ref¬ 
erences  22-41,  22-42)  having  varying  degrees  of  complexity,  with 
some  similarity  to  the  above  codes.  They  carry  many  species  and 
may  be  considered  nvltispecies  codes.  Table  22-1  contains  some 
important  E  and  F  region  species  which  are  carried  by  multispecies 
codes,  especially  by  the  NRL  Master  Code.  Multispecies  codes 
which  calculate  the  electron  temperature  also  consider  additional 
excited  states,  which  are  short-lived  and  decay  radiative ly  but  are 
not  carried  individually.  They  are  accounted  for  in  the  electron- 
temperature  equation  as  cooling  processes.  Examples  of  these 
states  are: 

N2(D2£),  N2(AZD),  Nz(c3n),  N2(B3H),  etc. 
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The  significance  of  these  and  the  states  shown  in  Table  22-1  is  ob¬ 
vious  for  the  optical  and  infrared  calculations. 


Table  22-1 .  Some  species  of  importance  fa  E  and  F 
region  multiple  species  computer  codes. 


Neutrals* 

Neutral 
Excited  States 

Positive  Ions 

Positively  Charged 
Excited  States 

N2 

n2(a3<) 

N2+ 

0 

oOd) 

o(’$) 

o+ 

o 

CM  .CM 

O  O 

CM 

O 

C„(bV) 
*'  g' 

► 

+  CM 

o 

N 

n(2d) 

Npp) 

N+ 

n+('d) 

N+('s) 

NO 

+ 

O 

z 

Me** 

MeO** 

Me+ 

MeO+ 

NOTES: 

One  may  include  minor  species  such  as  COj,  Ar,  etc. 

**Me  denotes  a  metallic  atom;  MeO  denotes  a  metallic  oxide. 
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Some  results  of  deionisation  calculations  are  shown  in  Figures 
2.1-2  through  22- 5.  These  results  are  taken  from  the  calculations  of 
Lewis  and  Ali  (Reference  22-43)  who  have  utilized  the  NRL  Master 
Code.  In  these  calculations  the  nitrogen- vibrational  temperature 
(Reference  22-40)  is  computed  assuming  a  harmonic  oscillator  model 
(Reference  22-44).  However,  one  may  calculate  (References  22-38 
through  22-40)  the  time  histories  of  the  individual  vibrational  levels 
instead  of  a  vibrational  temperature.  Of  course,  this  adds  to  the 
complexity  of  the  code.  The  reason  for  carrying  individual  vibra¬ 
tional  states  lies  in  the  fact  that  each  vibrational  state  of  N£  has  a 
different  reaction  rate  for  the  formation  of  NO^  (Reaction  22-6),  For 
example,  the  v  -  4  state  has  a  rate  coefficient  which  is  at  least  two 
orders  of  magnitude  larger  (Reference  22-29)  than  the  v  =  1  or  v  =  0 
state.  Figure  22-6  shows  a  deionization  calculation  where  the  time 
histories  of  the  vibrational  levels  are  indicated. 


Some  remarks  concerning  the  results  in  Figures  22-2  through 
22-5  are  in  order:  Figures  22-2  and  22-3  show  the  deionization  cal¬ 
culations  for  a  given  initial  condition  at  an  altitude  of  200  km  where 
densities  of  N?,  O,  and  O2  are  3.25  *  10 ^ ,  3.216  *  10^,  and 
3.5  y  10^  cm"  ,  respectively.  No  molecular  iors  are  present  be¬ 
cause  of  the  complete  ionization  of  the  region.  Tl.e  electron  density 
decreases  by  a  factor  of  2  over  —300  seconds  due  to  radiative  re¬ 
combination  resulting  in  an  increase  in  the  N  and  O  concentrations. 

An  interesting  picture,  however,  appears  for  a  lower  degree  of 
ionization  at  the  same  altitude  (Figures  22-4  and  22-5).  In  this 
case  the  degree  of  ionization  is  —  1  2  percent.  N,  and  are  depleted 
within  -10  seconds  (Figure  22-4).  NO*  is  abundant  and  is  mainly 
responsible  for  the  removal  of  the  electrons  for  times  beyond  a  few 
seconds.  The  electron  density  decreases  by  —3  orders  of  magnitude 
over  a  period  of  300  seconds  because  of  the  presence  of  molecular 
ions  (cf,  the  initial  case).  The  increase  in  the  neutral  particle 
densities  arises  from  dissociative  recomhination.  The  early  time 
increase  in  the  metastable  neutral  particle  concentrations  is  due  to 
direct  electron  impact  excitation  and  the  dissociative  recombinations. 

It  is  obvious  that  such  a  multispecies  treatment  cannot,  at  least 
at  this  time,  be  carried  !.n  three-dimensional  magnetohydrodynamic 
codes.  However,  if  the  purpose  of  the  deionization  calculations  is 
determination  of  the  electron  density,  then  one  may  use  the  Simple 
Chemistry  Code  (Reference  22-40)  and  eliminate  most  of  the  meta- 
stable  states  by  accounting  realistically  for  their  effects  on  the 
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.  Electron  and  atomic-ion  densities  for  a  given  initial 
condition  at  200  km  altitude. 
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Figure  22-3.  Other  species  densities  for  the  same  point  as  in  Figure  22-2. 
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deionization  process.  Figures  22-7  and  22-8  show  a  comparison  be¬ 
tween  the  electron  density  as  calculated  by  the  Simple  Code  (smooth 
curves)  and  the  NRL  Master  Code  (XXX  and  OOO)  for  six  different 
ionization  cases.  As  an  example,  Figure  22-9  shows  a  comparison 
between  some  ionic  species  densities,  for  a  given  case  (Case  I  of 
Figure  22-7)  as  calculated  by  both  codes. 

From  these  examples  one  may  infer  that  simplification  in  the 
chemistry  is  possible.  However,  if  one  is  interested  in  -ne  optical 
and  infrared  emission,  simplification  can  hardly  be  justified. 

22.3.3  The  Disturbed  D-Region  Multispectes  Codes 

The  D-region  chemistry  is  more  complex  than  that  of  the  E  and 
F  regions,  because: 

(a)  The  D  region  is  denser,  which  enhances  the  occurrence  of 
three-boc';'  orocesses,  especially  the  three-body  attachment 
of  electrons  to  02  . 

(b)  The  D  region  has  many  more  minor  species  (see  Table  22-2) 
which  play  an  important  role  in  ion-cluster  formation. 

The  deionization  processes  discussed  in  Section  22.  3.1,  e.  g.  , 
dissociative  recombination,  ior-mdecule  interaction,  and  charge 
exchange,  which  are  responsible  for  the  removal  of  electrons  in  the 
E  and  F  regions  ara  equally  valid  in  the  D  region.  However,  several 
different  physical  processes  arise  in  the  D  region.  The  three-body 
attachment  of  the  electron  to  according  to  the  reaction: 

e  i  O  *  M  -  O"  4  M  (22-11) 

£  L. 

becomes  an  important  electron  removal  process  and  has  a  rate  coef¬ 
ficient  (Reference  22-15)  of  -TO--^  Cm^  sec-'  at  D-region  tempera¬ 
ture.  Reaction  (22-11),  the  charge  exchange: 

°2  +  °3  "*  °3  +  °2  ’  (22-12) 

and  the  dissociative  attachment  to  start  the  negative-ion  reaction 
sequence.  Figure  22-10  shows  a  schematic  diagram  (Reference 
22-46)  for  the  formation  of  negative-ion  clusters  through  many  pos¬ 
sible  reactions  of  the  negative  ions  with  H20,  C02,  and  NO^. 
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Figure  22-9.  Comparison  of  ion  density  calculations  for  Case  I  of  Figure  22-7 
(smooth  curves,  NRL,  Simple  Code;  XX  00,  NRL  Master  Code). 
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Figure  22-10.  Schemotic  representation  of  the  formation  of  negative  ions  in  ionized  oir.  The  solid  lines 

indicate  probable  reactions,  on  the  basis  of  experimental  measurements  of  the  rate  coefficients. 
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Table  22-2.  Some  species  of  importonce  to  the  Sub-D-  ond 
D-regions  multispecies  computer  codes. 


Ncutroll 

Neutroi 

Excited  States 

Negatives 

Positives 

ch4 

N(20' 

co; 

H'(H20)n  1-5 

CO 

hjaV) 

co;(h20} 

h4(h2o)(hoi 

co2 

O(’d) 

C°4* 

HfHjOHNji 

H 

O(’s) 

C04(H20) 

n‘ 

HNOj 

0*(oV  ) 

/  9 

NOj 

no" 

HN°3 

OJb't4) 

2  9 

no‘(h2o) 

NO‘(C02> 

HO 

n°3  * 

N°‘(H2°)n  N3 

r*4 

o 

X 

"WU-S 

NOfNj) 

H 

2 

0* 

NOj 

H2° 

°2 

S 

H2°2 

OjfHjO) 

O* 

N 

°3 

4 

°2 

NO 

°; 

o2(H2oi 

c 

O 

z 

oo  no’  * 

4 

°4 

f-i 

O 

z 

6 

N2 

NjO 

N2°5 

0 

°2 

°3 

. 

£Tht  difference  between  the  designations  OONO*  end  NC^  ii  explained  on  page  18A-7 
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The  formation  of  positive -ion  clusters,  on  the  other  hand,  starts 

with  the  formation  of  0+  according  to: 

4 

o!,  <  O  4  M  -  O'  +  M  ,  (22-13) 

1*  Li  *X 

which  in  turn  reacts  with  H  O; 

2 

04  +  H20  -  D*(H  O)  +  02  ,  (22-14) 

to  produce  the  singly  hydrated  molecular  oxygen  ion. 

A  detailed  scheme  shown  in  Figure  22-1’  indicates  the  possible 
mechanisms  for  he  formation  of  positive-ion  clusters  (Reference 
22-46), 


22.3.4  Exomples  of  D-Region  Deionfzotion  Colculotions 

Niles  (References  22-46  and  22-47),  and  Scheibe  (Reference 
22-48)  have  developed  disturbed  D-region  deionization  codes  which 
are  multispecies  with  hundreds  of  reactions.  These  codes  have 
been  used  to  carry  out  several  deionization  calculations  (References 
22-46  through  22-51)  for  the  beta  patch  and  cool  D  region  outside 
the  fireball.  These  calculations  have  not  included  transport.  Tr.is 
was  the  case  also  for  the  E-  and  F-i'ugion  calculations  presented 
earlier.  However,  codes  are  being  developed  which  include  trans¬ 
port  and  a  limited  number  of  reactions  (Reference  22-52).  Since 
for  most  of  the  species  in  the  D  region  the  time  constant  for  re¬ 
action  is  much  less  than  the  time  constant  for  diffusion,  the  omis  - 
sion  of  transport  should  not  greatly  affect  atmospheric  deionization. 
However,  this  is  not  true  for  the  E  and  F  regions. 

The  calculations  of  the  multispccies  codes  presented  here  are  for 
the  60-km  region,  where  UHF  attenuation  is  greatest  (Reference 
22-46),  These  calculations  utilize  the  model  atmosphere  given  in 
Table  22-3  and  the  analytic  expressions  in  Table  22-4  for  the  elec¬ 
tron  productions  produced  by  beta- ray  ionization. 

A  composite  of  the  electron-density  decays  (Reference  22-46)  in 
the  beta  patch  at  an  altitude  of  150  km  is  given  in  Figure  22-12. 

The  decay  curves  are  identified  by  the  altitude  (in  kilometers)  to 
which  the  curve  applies.  The  form  of  the  electron-density  decay 
divides  readily  into  two  temporal  regions.  The  electrons  created 
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Table  22-4.  Analytical  expressions  for  electron  production 
[q(e,  t)  =  f(t)  crts"^  sec"M. 


Altitude 

(km) 

f(t) 

(t  in  sec) 

80 

4.7  x  107(1  +  tf1,2(l  +  15e'f/2°) 

75 

9.3  x  107(1  +t)“K2(1  +12e“f/20) 

70 

1.4  x  108(1  +tfL2(l  9  e“f/2°) 

65 

1.6  x  108(1  +  t)"1,2(l  +  12  e_f/2°) 

60 

1.8  x  108(1  +  tf1,2(l  +  15e‘t/2°) 

55 

1.6  x  108(1  +t)“1,2(l  +  12e'f/2°) 

50 

1.4  x  108(1  +  t)“1*2(l  +9e'f/2°) 

by  prompt  radiation  decay  initially  by  dissociative  recombination 
with  N2  and  then  with  C>2  in  Region  1.  Region  II  begins  when  the 
electron  density  is  affected  by  the  ionization  produced  by  the  beta 
rayu.  In  Region  II,  the  electron  density  is  lost  primarily  by  dis¬ 
sociative  recombination  with  NO+  in  the  upper  part  of  the  beta 
patch  and  by  electron  attachmenf  to  in  the  lower  part.  In  agree¬ 
ment  with  this,  the  dependence  of  the  electron  density  on  the 
source  term,  q(e,t),  changes  from  square  root  in  the  upper  part  to 
linear  in  the  lower  part.  In  between,  the  electron  density  decay  is 
complex. 

In  Region  1,  the  decay  of  electron  density  is  given  by: 

dn  /dt  =  Oc  n  n  ,  (22-15) 

e  d  e  + 

where  n+  is  the  total  positive-ion  density  and  0^  is  the  lumped  dis- 
sociative -recombination  coefficient.  Since  ne  *•-  n+  ,  the  solution 
can  be  written  as: 
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Figure  22-12.  Decoy  of  electron  density  in  the  beta  patch  at  150  km 
altitude.  The  numbers  associated  with  the  decay  curves 
give  the  altitude  in  kilometers  (Reference  22-46). 
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n  =  (n  )  /[(l  +  (n  ,  (22-16) 

e  e  0  e  0  d 

where  (n  )  is  the  initial  electron  density  immediately  after  the 
burst.  6  ° 

Since  the  dissociative-recombination  rate  constants  for  Nt  and  Ot 

*7  o  |  c,  l 

recombining  with  electrons  are  about  3  x  10"'  sec-1  and  the 

initial  electron  density  is  about  10^  cm"^,  Equation  (22-16)  reduces, 
at  t  >  10"^  sec,  to: 

n  =  1/at  .  (22-17) 

e  d 

This  says  that  if  there  were  no  source  term  and  if  chemistry  were 
unimportant,  then  the  electron-density  decay  would  be  the  same  for 
all  altitudes  and  all  levels  of  initial  ionization  in  the  beta  patch. 

Figure  22-12  shows  a  temporal  region  between  10"^  and  10"^  sec 
which  corresponds  to  this  situation,  within  which  the  time  required 
for  the  electron  density  to  decay  to  a  specified  value  apparently 
could  be  readily  determined.  However,  as  seen  in  Figure  22-12, 
this  is  not  the  case. 

In  Region  II,  the  decay  of  electron  density  in  the  upper  part  of  the 
beta  patch  is  given  by: 

dn  /' dt  =  q(e,  t)  -  Q  n  n  .  (22-18) 

e  d  e  ‘t 

The  quasi-equilibrium  solution  when  n  ~  n4  ^s: 

n  =  |q(e,t)/a]0,5  .  (22-19) 

e  a 

In  the  lower  part  of  the  beta  patch  the  decay  of  electron  density  is 
given  by: 

dn  /dt  -  q(e,t)  -  An  ,  (22-20) 

e  e 

where  A  is  the  electron  attachment  frequency.  Since  A  is  deter¬ 
mined  primarily  by  three-body  attachment  of  O 2  and  since  the  O2 
density  is  nearly  time-invariant,  the  quasi-equilibrium  solution  is: 

n  -  q(e,t)/A  .  (22-21) 

e 
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The  relative  importance  of  the  various  reactions  which  affect  the 
electron  density  changes  with  time.  The  rates  for  the  populating  re¬ 
actions  at  60  km  are  shown  in  Figure  22-13,  where  Reactions  1,  15, 
and  201  denote,  respectively,  photodetachment,  associative  detach¬ 
ment,  and  collisional  de-excitation  detachment  reactions  of  . 

The  depopulating  reactions  for  the  electron  density  giving  rise  to 
atmospheric  deionization  are  shown  in  Figure  22-14  along  with  their 
rates.  It  is  seen  that  dissociative  recombination,  first  with  N*  and 
then  with  ,  is  the  dominant  loss  process  out  to  about  10  sec,  at 
which  time  three -body  attachment  to  0£  becomes  dominant. 

The  importance  of  cluster  ions  can  be  deduced  from  the  rates  for 
Reactions  296,  327, and  334  (Figure  22-14).  These  represent,  respec¬ 
tively,  the  dissociative  recombinations  of  O^,  H^O+ (H20),  and  H^O^HjOlj 
with  electrons.  Since  for  this  altitude  the  initial  conditions  for 
the  calculations  assumed  the  ambient  densities  of  these  ions  to  be  zero, 
the  increase  in  rates  reflects  the  rapid  growth  in  cluster-ion  densities. 

Figures  22-15  through  22-17  show  the  results  (Reference  22-46) 
of  a  deionization  calculation  for  a  beta  patch  at  60  km.  The  figures 
show  the  negative  ions,  the  positive  ions,  and  the  neutral  densities, 
respectively, 

Scheibe  (Reference  22-53)  has  produced  a  simplified  analytic 
model  for  the  disturbed  D-region  chemistry  by  obtaining  analytical 
approximations  for  the  species  densities.  His  approach  is  based  on 
the  fact  that  for  a  moderate  to  high-level  ionization  without  appreci¬ 
able  heating,  a  sudden  ionization  phenomenon  can  be  separated  into 
four  distinct  regions  in  time.  The  first  period  comprises  the 
ionization,  dissociation,  and  excitation  of  the  air  by  x-rays,  and  is 
over  in  a  very  short  time.  The  second  phase  includes  recombination, 
charge-exchange,  and  ion-neutral  reactions  producing  NO,  O,  and 
N.  Even  though  the  time  scale  for  phase  two  is  long  compared  to 
phase  one,  it  is  usually  complete  in~l  sec.  At  this  stage  the 
charged-particle  chemistry  is  no  longer  coupled  to  the  neutrals 
chemistry.  The  third  phase  involves  the  decay  of  O,  N,  02<a  lA), 
and  NO  and  the  production  of  O^  and  N02.  The  final  phase  is  the 
late-time  deionization  whose  solution  is  based  on  the  fact  that 
after  a  short  time  the  charged-pai tide  concentrations  will  achieve 
a  steady-state  condition  and  the  concentrations  can  be  obtained  by 
algebraic  relations.  However,  this  may  not  always  be  true  for  such 
species  as  NO\  NO",  CO",  and  e.  For  these  species  one  could 
solve  ihe  appropriate  rate  equations. 
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Figure  22-13.  Electron  populating  rates  ot  60  km  altitude  in  the  beta 
potch  (Reference  22-46). 
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Figure  22-14.  Electron  depopulating  rotes  at  60  km  altitude  in  the  beta 
patch  (Reference  22-46). 
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Figure  22-15.  Electron  and  negative-ion  densities  at  60  km  altitude  in 
the  beta  patch  (Reference  22-46). 
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Figures  22-18  through  22-20  show  the  comparison  of  certain 
species  densities,  for  a  high-level  ionization,  calculated  by  a  de¬ 
tailed  rrultispecies  code  (solid  lines)  and  analytically  (Reference 
22-53)  (□□,  CO,  A).  Figures  22-21  through  22-23  show  the  same 
comparison  for  a  low-level  ionization. 

The  results  show  that  the  analytic  model  is  fairly  good  for  a  high 
ionization  level  but  not  as  good  for  low  ionization  levels. 

22.4  LUMPED-PARAMETER  METHOD 

The  utilization  of  multispecies  codes  can  involve  prohibitive 
levels  of  effort,  at  least  at  the  present  stage  of  code  and  machine  de¬ 
velopment,  in  problems  where  deionization  calculations  are  re¬ 
quired  for  a  very  large  number  of  points  in  space  and  time.  There¬ 
fore,  approximations  are  in  order  for  a  large  group  of  problems 
where  one  can  thereby  reduce  computational  time  requirements  ap¬ 
preciably.  Several  approaches  were  discussed  earlier,  e.g.,  the 
elimination  of  certain  species  and  reactions  which  have  negligible 
effects  on  the  quantities  being  computed.  For  electron  density  calcu¬ 
lations,  Lewis  and  Ali  (Reference  22-43)  have  shown  that  for  the  E 
and  F  regions,  one  can  obtain  accurate  results  (less  than  20  per¬ 
cent  error)  using  a  simpler  code  (See  Figures  22-7  and  22-8).  For 
the  D  region,  Scheibe  (Reference  22-53)  has  decoupled  the  charged 
chemistry  from  the  neutral-species  chemistry  and  has  used  steady- 
state  relations  to  obtain  the  appropriate  concentrations.  For 
electron-density  calculations  his  approximations  give  rather  good 
results,  at  least  for  a  high  level  of  ionization  (Figure  22-20). 
Simplified  D- region  deionization  models  for  nuclear  environments 
based  on  eliminating  specified  species  and  reactions  have  also  been 
described  by  Bortner  and  Galbraith  (Reference  22-49),  and  Ory 
and  Gilmore  (Reference  22-54). 

Another  approach  to  the  reduction  of  deionization  computation 
times  is  to  lump  those  species  together  for  which  effective  rate  co¬ 
efficients  can  be  determined.  This  method  has  been  used  extensively 
in  systems  studies.  Knapp  (Reference  22-55)  has  recently  developed 
a  model  for  the  disturbed  D  region  using  the  lumped-parameter 
method. 

Deionization  models  in  which  species  are  lumped  together  and 
effective  reaction  rate  constants  are  used  have  been  described  ex¬ 
tensively  in  the  literature  (References  22-56  through  22-62).  For 
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Figure  22-18.  Sixty-km,  high-level  ionization  (imooth  curve,  detoited 
species  code; CO  Au,  analytic  model)  (Reference  22-53' 
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Figure  22-19.  Sixty-km,  high-level  ionizotion  (smooth  curve,  detailed  multi¬ 
species  code;  0  □  ,  onolytic  model)  (Reference  22-53). 
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Figure  22-21.  Sixty-km,  low-level  ionizotion  (smooth  curve,  detoiled  mult! 

species  code;  ADO,  onolytic  model)  (Reference  22-53). 
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Figure  22-22.  Sixty-km,  low-level  ionization  (smooth  curve,  detailed  multi¬ 
species  code;  C  □,  analytic  model)  (Reference  22-53). 


Rtvitton  No,  2,  M«v  1974 


22-39 


DENSITY  {cm" 


CHAPTER  22 


altitudes  below  about  100  km  a  three-charged-species  model  has 
been  used  for  radar  and  communication  system  studies.  In  this 
model  a  single  representative  positive-ion  species  and  a  single 
representative  negative-ion  species  are  used  and  the  neutral  species 
are  considered  to  comprise  an  unchanging  background.  The  reac¬ 
tions  considered  are  electron  attachment,  electron  detachment, 
electron  recombination  to  positive  ions,  and  ion-ion  mutual  neutrali¬ 
zation.  The  differential  equations  describing  the  species  densities 
are: 


dn 

—7 —  =  q  -  Of  n  n  -  An  +  Dn  (22-22) 

dt  d  e  •*  e 

dn 

dt  i  -  +  e  ' 

dn 

— —  -  q  -  Of  n  n  -  a.n  n  (22-24) 

dt  4  d  e  t  i  -  + 

where ; 

n  =  number  density  of  electrons 
e 

n  =  number  density  of  positive  ions 

n  =  number  density  of  negative  ions 

A  -  attachment  rate  of  electrons  to  neutrals 

D  -  detachment  rate  of  electrons  from  negative  iens 

0^  -  recombination  rate  coefficient  of  electrons  and 

positive  ions 

a.  mutual  neutralization  rate  coefficient  of  positive 
and  negative  ions 

q  =  electron-positive  ion  production  rate. 

All  reaction  rates  are  functions  of  temperature  and  neutral  species 
densities  at  least.  The  ionization  sou.ee  q  is  a  function  of  time  and 
space.  Equations  (22-22)  through  (22-24)  represent  only  two  inde¬ 
pendent  relations  since  charge  neutrality  (n  =  n  *  n  )  must  be 
preserved. 
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Since  the  lumped  rate  coefficients  comprise  a  combination  of  in¬ 
dividual  rate  coefficients  for  the  large  number  of  species  actually 
present,  the  general  problem  of  finding  overall  reaction  rates  as  a 
function  of  time  requires  solution  of  the  complete  set  of  differential 
equations.  The  lumped -parameter  model  is  most  useful  when  the 
physical  approximation  is  made  that  the  lumped  rate  coefficients  are 
independent  of  time  and  initial  conditions.  This  approximation  can 
only  be  justified  either  by  determining  from  detailed  multispecies 
calculations  that  for  a  specified  range  of  conditions  the  lumped  co¬ 
efficients  are  essentially  constant  or  by  determining  a  set  of  coef¬ 
ficients  that  provide  empirical  agreement  with  atmospheric  data  for  a 
specified  range  of  conditions.  The  latter  approach  has  been  used  in 
determining  lumped-parameter  reaction-rate  models  where  elec¬ 
tron  density  and  total  ion  density  were  the  principal  quantities  of 
interest. 

While  the  lumped-parameter  model  has  significant  limitations  and 
must  be  used  with  caution,  its  use  can  be  justified  in  a  number  of 
applications.  It  has  been  possible  to  derive  reaction  rates  which  pro¬ 
vide  reasonable  agreement  with  detailed  calculations  and  ionization 
conditions  of  interest  in  radar  and  communication  system  studies. 

The  lumped-parameter  model  is  simple  enough  to  provide  a 
physical  understanding  of  the  ionization  as  a  function  of  altitude 
and  ionization  conditions.  Mathematical  approximations  can  be 
made  which  provide  solutions  to  the  lumped-parameter  differential 
equations  which  are  in  some  cases  thousands  or  even  millions 
of  times  faster  than  the  complete  solution  of  the  multispecies 
equations . 

The  principal  limitation  of  the  lumped-parameter  model  is  its 
range  of  applicability.  When  the  reaction  rates  are  assumed  to  be 
only  a  function  of  altitude  the  model  v  ill  not  account  for  changes  in 
the  dominant  ion  species  with  ionization  conditions.  The  formation 
of  minor  neutral  species  during  the  initial  energy  deposition  and  sub¬ 
sequent  deionization  can  significantly  modify  the  effective  recombina¬ 
tion  and  detachment  rates. 

The  three-ion-species  differential  equations  can  be  solved  in 
closed  form  if  it  is  assumed  that  ft  -  ft^  =  ft.  and  q  is  a  constant. 

For  the  general  case  where  ftj  t  ft[  approximate  mathematical  solu¬ 
tions  of  the  differential  equations  have  been  developed  by  Bogusch 
(Reference  22-59)  and  by  Ek  (Reference  22-i'»0).  For  conditions 
where  q  varies  slowly  enough  with  time  so  that  after  the  initial 
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transients  die  out  the  ion  densities  approach  a  condition  where  they 
are  essentially  in  equilibrium  with  a  at  any  given  instant  the  follow¬ 
ing  quasi-equilibrium  solutions  can  be  used: 


(22-25) 


(22-26) 


(22-27) 


The  relation  for  determining  the  effective  recombination  coefficient 
<*e (a  cubic  in  ft  )  may  be  solved  directly,  or  an  iterative  procedure 
may  be  used.  The  quasi-equilibrium  solution  is  often  written  in 
terms  of  X,  the  ratio  of  negative-ion  to  electron  densities  instead  of 
fte  (Referen<e  22-63): 


n  -  (1  ♦  A)n 

e 


n 


(22-28) 


(22-2  0) 


X 


_ A 

D  ‘  0.(1  *  X)n 
i  e 


(22-30) 


Table  22-5  shows  reaction  rates  and  coefficients  chosen  for  use  with 
the  lumped-parameter  model  in  appropriate  studies  (Reference 
22-61).  The  values  giver,  are  consistent  with  laboratory  measure¬ 
ments  for  diatomic  ions.  As  discussed  in  Chapters  2  and  5  the 
dominant  atmospheric  ions  below  about  80  km  for  moderate  to  low 
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ionization  levels  are  hydrated  ions.  Since  hydrated  ions  have  much 
larger  rate  constants  for  dissociative  recombination,  it  would  be  ex¬ 
pected  that  the  reaction  rates  given  in  Table  22-5  would  only  be  ap¬ 
plicable  for  moderate  to  high  ionization  levels.  Figure  22-24  shows 
comparisons  of  electron  densities  obtained  using  the  lumped- 
parameter  model  with  those  obtained  through  a  multispecies  code 
A  comparison  of  the  one-way  vertical  absorption  resulting  from 
these  electron  densities  for  UliF  propagation  is  shown  in  Table 
22-6.  As  expected,  the  lumped-parameter  model  gives  good  re¬ 
sults  for  the  first  few  hundred  seconds  and  then  overestimates 
electron  densities,  owing  to  the  formation  of  hydrated  ions.  The 
comparison  shown  is  not  intended  to  imply  that  similar  agreement 
will  occur  for  other  ionization  conditions,  since  a  number  of  com¬ 
pensating  factors  are  involved  which  depend  on  the  level  of  ioniza¬ 
tion  and  the  altitude  of  interest. 


Table  22-6.  One-way  vertical  attenuation  in  dB  for  UHF  radars. 


1 

10 

32 

100 

316 

1 000  sec 

DETAILED  CODE  221 

50 

11 

D 

0.94 

0.35 

SIMPLE  CODE  271 

48 

11 

B 

1.0 

0.47 

A  luniped-parameter  deionization  model  in  which  the  reaction 
rates  are  dependent  on  ionization  conditions  has  been  developed  for 
use  in  communication  system  computer  codes.  The  model,  ob¬ 
tained  by  fitting  data  obtained  from  multispecies  calculations,  is 
primarily  applicable  for  cases  where  the  minor  neutral  species  com¬ 
position  is  not  significantly  modified  by  deionization  processes,  i.e.  , 
initial  ionization  less  than  about  10^  ion-pairs  cm-3  and  ion-pair  pro¬ 
duction  rates  less  than  about  10^  ion-pairs  cm"^  sec"^.  Lumped- 
parameter  models  for  quasi-equilibrium  conditions  which  include 
several  positive  and  negative  ions  and  effective  reaction  rates  which 
are  functions  of  both  the  ion-pair  production  rate  and  the  dominant 
minor  neutral  species  are  being  studied. 
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Figure  22-24.  Beta-patch  electron-density  profiles  at  stated  times  as 

calculated  by  detailed  or  simple  cades  (Reference  22-46). 

The  relationship  between  the  ion-pair  production  rate  and  the 
quas i- equilibrium  electron  density  in  the  D  region  is  often  written  in 
the  form: 


where  :  is  knoun  various!)  as  the  steady-state  coefficient,  effective 
recombination  coefficient,  or  effective  electron  loss  rate.  This  form 
is  particularly  useful  in  comparing  atmospheric  measurements  where 
determinations  of  q  and  nc  are  made.  The  quantity  '  can  be  ex¬ 
pressed  in  terms  of  X  or  0!(,  as  follows: 

(0!  *  Xtt. Ml  •  XI  ;  (22-321 

d  i 
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Figure  22- 25  shows  calculations  of  Y  obtained  from  detailed 
multispecies  code  data  for  parametric  values  of  q.  At  night  v  has 
a  relatively  strong  dependence  on  q.  Also  shown  are  atmospheric 
measurements  reported  by  Megill  et  al  (Reference  22-621  for  day¬ 
time  polar-cap  absorption  (PCA)*  events  (q  between  10^  and  103  ion- 
pairs  cm"3  sec"l).  The  differences  in  the  theoretical  and  measured 
values  of  •  are  believed  primarily  due  to  missing  reactions  in  the 
positive-ion  cluster  chain  starting  with  NO4,  and  to  differences 
between  the  minor  neutral  species  concentrations  used  in  the  multi¬ 
species  code  calculations  and  those  actually  present  at  high  altitudes. 

For  altitudes  above  about  100  km  both  atomic  and  molecular  posi¬ 
tive  ions  are  important  and  must  be  modeled  separately  owing  to 
significantly  different  Loss  rates.  Simplified  models  which  use  a 
representative  molecular  positive  ion  and  either  a  representative 
atomic  positive  ion  or  explicit  models  of  O4  and  N4  have  been  de¬ 
veloped  for  use  in  radar  systems  computer  codes.  The  reaction 
rates  are  formulated  in  terms  of  a  molecular  nitrogen  vibrational 
temperature  and  an  electron  temperature  which  are  determined  from 
approximate  models  for  the  natural  atmosphere  and  the  amount  of 
energy  deposited  during  ioni/.ation.  The  problem  of  determining  de¬ 
ionization  in  the  E  and  F  regions  is  complicated  by  atmospheric 
motion  and  the  effect  of  deionization  on  the  atomic  and  molecular 


3  •! 

t  <tm  IK 

Figure  22-25.  Multispecies  code  colculotions  of  . 

••'Quite  recently  this  phenomenon  has  come  to  be  known  as  the  Solar 
Proton  Event  (SPE) . 
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neutral-species  densities.  Generally,  even  when  a  lumped- 
parameter  model  is  used,  the  deionization  must  be  solved  sequen¬ 
tially  for  a  series  of  time  steps.  Nevertheless,  for  certain  applica¬ 
tions  the  lumped-parameter  models  are  reasonably  accurate  and 
are  significantly  faster  than  numerical  solutions  of  the  differential 
equations . 
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